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Abstract 

The capacitance of arrays of vertical wrapped-gate InAs nanowires are analyzed. With the 
help of a Poisson-Schrodinger solver, information about the doping density can be obtained di- 
rectly. Further features in the measured capacitance-voltage characteristics can be attributed to 
the presence of surface states as well as the coexistence of electrons and holes in the wire. For 
both scenarios, quantitative estimates are provided. It is furthermore shown that the difference 
between the actual capacitance and the geometrical limit is quite large, and depends strongly on 
the nanowire material. 
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I. INTRODUCTION 

Large efforts are currently made to develop new schemes for the manufacturing of 
nanowire (NW) devices fl. E.g., recent progress in vertically processed semiconductor NWs 
with a cylindrically symmetric gate electrode, such as the wrapped insulator-gate field effect 
transistor have shown great promise for application in future highly scaled electronic 
devices. Besides the possibility of incorporating wrapped gates for improved gate control, 
NWs offer other intriguing possibilities not easily accessible to planar designs, such as an 
inherent one-dimensionality as well as relaxed constraints in terms of combining highly 
lattice-mismatched materials along the NW channel due to radial relaxation of interface 



strain 



In order to facilitate the development of high performance NW electronic devices, it 
is important to accurately characterize the doping profile and interface properties of the 
system as well as further fundamental material parameters. Si-based MOS-stacks have long 
benefited from well developed capacitance-voltage (CV) measurement schemes for device 
characterization. However, similar methods have so far been largely unavailable for NW 
devices. Recently CV characteristics of individual Ge NWs on SiO^/Si substrates were 



measured 
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Here we are focusing on vertical wrapped gate NWs for which we have earlier 
demonstrated a scalable processing protocol for routine CV spectroscopy [6]. An overview of 
the device structure considered here is illustrated in Fig. [IJ The device consists of an 11x11 
array of InAs NWs, covered with a 10 nm conformal Hf02 dielectric layer and a top Cr/Au 
metalization as one of the capacitor electrodes. The capacitance is measured between this 
gate electrode and the core of the wire, which is connected via the substrate. Details on the 
fabrication protocol and measurement techniques can be found elsewhere js]. In this paper we 
further develop the quantitative as well as qualitative understanding of these measurements, 
and also discuss various features of the CV-characteristics in detail. 



II. THE MODEL 

The nanowires are modelled as cylinders where the radius of the InAs will be referred to 
as -R nw , whereas -R ga t e is the radius of the InAs and Hf0 2 together. Typical measures are 
length h = 680 nm, and radial dimensions R QW = 27 nm, -R ga t c = 37 nm. 
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FIG. 1: Illustration of the array of NWs. Each NW consists of an InAs core surrounded by a 
dielectric layer of Hf0"2 and a Cr/Au gate. The lifting layer of polymer is inserted to reduce the 
capacitance between gate and substrate. 



We assume that the unintentional doping of the NWs is of n-type. By changing the 
applied gate voltage, V^ ate , the charge within the NW, Q, can be affected. This defines the 
capacitance 



« ^gate 

The total charge Q is the spatial integral over the charge distribution p(r) containing 
three components p = e{N£> — n e + rih) inside the InAs nanowire. Here, e is the positive 
elementary charge and Njj, n e and are the concentrations of ionized donors, electrons, 
and holes, respectively. The donors are assumed to be entirely ionized at room temperature 
as considered throughout this paper. All concentrations are assumed to depend only on 
the radial distance r from the centre of the NW. This greatly facilitates the calculations 
as it reduces the problem to one dimension. Working with homogenous doping profiles is 
motivated by averaging over arrays of long wires. However, microscopic fluctuations are 
neglected by this approach. 

The effective density of states masses for electrons and holes in InAs are m* = 0.023m e 



and m* h = 0.41m e respectively [3, Qj- Non-parabolicity effects are neglected albeit they can 
cause some quantitative deviations for large positive biases. For the determination of hole 
concentrations, we use the bandgap E g = 0.54eV for InAs NWs with wurtzite structure |9[]. 
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ie relative dielectric constant is e ~ 15 both for the InAs NW and the HfOo insulator 
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To calculate the capacitance, a Poisson-Schrodinger solver was implemented similar to the 
one demonstrated in Q,Q. For a given charge distribution p(r), the cylindrical geometry 
permits the use of Gauss law to solve Poisson's equation for the electrostatic potential 0(r). 
Thus, we obtain the Hamiltonian 

H = -^—V 2 -e<P(r) (2) 
2 ml 

together with the boundary condition for the wavefunction ^/(R nw ) = assuming that the 
Hf02 layer constitutes an infinite barrier. Due to the cylindrical symmetry of the problem, 
Bessel functions were used as a basis set for the diagonalisation of Eq. ([2]), resulting in the 
wave functions ^ U j(r, tp) and corresponding energies E U j. Here v and j are the angular and 
radial quantum numbers, respectively. The relatively long length of the wires compared to 
their radial dimensions justifies the use of a ID density of states along the z-direction, g(E z ). 
Using Fermi-Dirac statistics the electron concentration is then given by 

n e (r) = £ |^(r)| 2 f° 1 (3) 

^ Jo exp((E z + E UJ )/kT) + 1 

Here the electrochemical potential Ep was set to zero in the wire, which is used as a reference 
point throughout this paper. Together with the similarly calculated hole concentration 
Uh, this provides a new charge distribution p(r). By iterative solution of the Poisson and 
Schrodinger equation, a self-consistent solution is obtained. 

In Fig. [2] calculated potentials and electron distributions are shown. Figs. [2^) and b) 
correspond to a positive electrical potential at the gate electrode. The corresponding bending 
of the conduction band edge causes electrons to accumulate close to the semiconductor- 
oxide interface (accumulation mode). Figs. [2b) and d) correspond to a negative electrical 
potential at the gate electrode. Here the electrons are depleted close to the surface of the 
wire (depletion mode). 

In the simulation we assume that the insulating oxide layer is neutral of charge and use 
the corresponding electric potential </>o(-Rgate) at the gate as a parameter. However, trapped 
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FIG. 2: a) Potential V(r) = —exf>{r) and valence band edge V(r) — E g for <^o(-Rgatc) = 1.08 V, 
corresponding to the operation point 1 of Fig. El The dashed line depicts Ep of the wire, which is 
set to zero as reference, b) corresponding electron distribution n e (r). c) and d) are the same plots 
for 4>o{Rg a te) = —0.46 V, (operation point 2). Parameters: Np> = 2.2 • 10 18 /cm 3 [dashed line in b)], 
R nw = 27 nm. 

charges are likely to occur both at the surface of the wire and within the oxidejjjj]. If an 
areal charge density a is present at position -R trap the potential for r > -R t ra P is modified 
according to 

(4) 



(r) = 0o( r ) ~ a ~^ tra P ] n 



££q \ J Hrap, 

while 4>{r) is not affected for r < i2t rap . Experimentally, one measures the applied gate 
voltage Vgate which is related to the actual electric potential at the gate via 

Wcr ~ XlnAsNW 



gate 



gate . 



+ 



(5) 



Here Wc T is the work function of Cr and XlnAsNW the electron affinity of the InAs in the 
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nanowire. As long as the trapped charges are unchanged they provide a constant bias offset 



^gate = </>o(-Rgate) + Vq (6) 

between the measured bias and the parameter 0o(-Rgate) used in the simulation. As both 
the magnitude of the trapped charges and the electron affinity of InAs in the wurtzite 
modification appearing in the nanowire are not a priori known, we use Vq as a fit parameter 
in the subsequent data analysis. 



III. UNDERSTANDING THE EXPERIMENTAL DATA 

In Fig. [3] the experimentally measured capacitance for the device with mean radius 27 nm 
is shown together with our calculations for different doping densities. The main trend is the 
increase of capacitance with bias. This can be attributed to the location of charges inside 
the NW. From Gauss law it follows that the change in Q induced by a change in bias, will 
be larger if p(r) is affected at large values of r. As can be seen in Fig. [2b) (corresponding to 
the operation point 1 in Fig. [3]), the charges are mainly located close to the NW surface in 
the accumulation region. In contrast, they are close to the centre of the wire at operation 
point 2, see Fig. [2U). This explains, why the capacitance is higher at point 1 than at point 
2. 

The measurements reported here are performed by adding a small AC-signal with a 
frequency of 20 MHz atop the DC gate-bias. The impedance is found to be almost com- 
pletely imaginary showing a negligible influence of any series resistance to the accumulation 
capacitance. For positive bias the measured capacitance shows less than 10 % frequency 
dispersion between 50 kHz and 100 MHz demonstrating good properties of the metal-oxide- 
semiconductor structure and that a reliable quantitative data analysis can be performed for 
V^ at e > 0. In contrast, for V^ ate < the frequency dispersion is more pronounced, which we 
attribute to the presence of surface states, see Sec. IIII CI 

The calculations provide a capacitance per wire length, as given in Fig. [31 For the exper- 
imentally measured capacitance C, the gate substrate capacitance, which is independent of 
bias, see p], has been subtracted and the result has been divided by an effective length L, 
which has been taken as a second fit parameter in addition to Vq. The data shown assumes 
an effective length corresponding to an effective number of 90 wires with the nominal length 
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FIG. 3: Experimentally measured capacitance and theoretical fits for a mean wire radius R QW = 27 
nm. The green curve shows the measured capacitance when sweeping from positive to negative 
bias, while the yellow curve is for the opposite sweep direction. The gate-substrate capacitance, 
which is independent of bias, has been subtracted. 

of 680 nm. This has to be compared with a total of 121 wires in the fabricated array. The 
difference can be attributed to lacking electrical connection of some wires and/or variations 
from the nominal length. 

Each fabricated array showed a certain distribution of radii -Rnwjf]]. To accommodate for 
this, an ensemble of three different values of R nw was used in the simulation for each device, 
while the oxide thickness of 10 nm was held constant. The calculated curves in Fig. [3] and H] 
are the results from such ensembles with appropriate weights for each radius reflecting the 



measured radius distribution. 



In Fig.[3]it is indicated that different factors, such as geometry, doping, surface states, and 
inversion (holes), affect the capacitance at different gate voltages. We will in turn explain 
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i? nw (nm) 


# wires 


N D = 10 18 /cm 3 


offset V (V) 
N D = 2.2 • 10 18 /cm 3 


Nd = 4 • 10 18 /cm 3 


25 


71 


0.27 


0.42 


0.63 


27 


90 


0.26 


0.39 


0.60 


28.5 


86 


0.13 


0.32 


0.53 


30 


90 


0.12 


0.28 


0.51 


31 


84 


-0.05 


0.14 


0.38 



TABLE I: Effective number of wires and bias offset V$ resulting in a best fit of the capacitance-bias 
characteristics for different arrays of NWs and different assumed doping densities. 

the contribution from each factor. Measurements indicate that the down-sweep curve (green 
curve) is closer to equilibrium as time-dependent sweeps are more stable and the up-sweep 
(yellow curve) relaxes towards the down-sweep curve. Therefore, all fits of our stationary 
model are made to the down-sweep curve. 

A. Geometry 

According to Fig. [3] the capacitance is independent of the donor concentration in the accu- 
mulation mode, but rather depends on the geometry of the wire. The reason is that during 
accumulation electrons are added close to the semiconductor-oxide interface, see Fig. [2Jd. 
Thus different doping concentrations result in the same capacitance. The capacitance in 
this region can therefore not be used to determine the doping concentration. Fitting to the 
experimental data instead determines the effective length of the wires. 

In Fig.HJthe capacitance (experimental down-sweep) for NWs with different radii R nw are 
shown. All fabricated arrays have 121 NWs, and the effective number of wires used for fitting 
L is displayed in Table H As a general trend, the capacitance increases with the radius R nw 
as there is a larger surface area of the NWs. Note that the capacitance in the accumulation 
region is significantly lower than the geometrical capacitance C/L = 2neeo/ ln(i?g ate /-Rnw), 
which is 2.48 ■ 10~ 18 F/nm for R nw = 25 nm and -R ga te = 35 nm. This deviation is strongly 
dependent on the density of states for the nanowire wire material, see Sec. IIVI 
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FIG. 4: Experimentally measured down-sweep capacitance (green) and simulation results for dif- 
ferent NW radii R nw . The doping densities for the different simulations are identical to Fig. 

B. Doping 

As the bias is lowered, the electrons will be depleted in the NW. A large doping con- 
centration provides a positive background charge with a strong attraction for the electrons, 
requiring a higher negative gate voltage to deplete all the electrons. Thus, the doping 
concentration strongly affects the slope of the capacitance curve around V^ at0 = V. Com- 
parison between fits and experimental data shows that Nq = 4 • 10 18 cm -3 is too flat while 
= 1 • 10 18 cm -3 is too steep. The doping of = 2.2 • 10 18 cm~ 3 fits very well up to 
the "kink" at point 2 of the experimental down-sweep curve in Fig. [3j This also holds for 
all devices shown in Fig. HI As all devices were manufactured using the same method, this 
provides a consistent estimate of 2.2 • 10 18 cm -3 for the doping density. 
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C. Surface States 



The experimental curves exhibit a "kink" at point 2 of Fig. [3l which is most clearly seen 
for the down-sweep curve. For biases below the kink the measured capacitance is significantly 
larger than the simulation result. In this section we attribute this to the emptying of surface 
states, situated at the semiconductor-oxide interface. 



The surface states are assumed to be located in the forbidden gap 15fl, and are treated 

n 

as donors [16(. As the states are located in the bandgap, they remain filled as long as Ep 
is in the conduction band. Theoretical evaluation of the potential at -R nw , for the doping 
of N D = 2.2 • 10 18 cm~ 3 , shows that E F is 0.1-0.2eV below the conduction band edge at the 
kink, see Fig. [2b). This makes it feasible that the kink is an effect of such surface states. 
For gate voltages more negative than the one corresponding to the kink, surface states are 
emptied which produces the discrepancy between simulations and experiment. 

While conduction band electrons as majority carriers can almost adiabatically follow the 
AC signal at 20 MHz, the filling and emptying of surface states is a slower process. Therefore 
these states do not directly contribute to the capacitance under our measurement conditions. 
For lower AC frequencies, we observe an enhancement of the capacitance in the bias region 
below the kink which we attribute to the direct contribution of surface states. The situation 
is however intricate due to a strong variation of charging times and its analysis is a matter 
of ongoing research. Thus we focus on the high-frequency data at 20 MHz here. 

Even if the surface states do not directly contribute to the capacitance signal, they affect 
the capacitance- voltage characteristics indirectly by the presence of a bias-dependent surface 
charge. This positive charge screens effectively the gate bias from the NW, so that a more 
negative V gSute is needed to deplete the electrons from the wire. This can be quantified 
using Eqs. (Hf5|) . An additional (bias-dependent) surface charge density 8a at R nm gives an 
additional shift of the gate bias V ga _ tc = 4>o(R ga ,t c ) + V + 5V with 

eeq \ R nw J 

This explains the shift of the capacitance curve to lower bias while emptying the surface 
states. With 5V being of the order of —1 V, positive surface charge densities of ~ 10 13 /cm 2 
are observed for gate biases far in the depletion region. 

To summarize, if surface states are not included in the model, the simulations should 
be fitted to the part of the down-sweep curve right of the kink. The difference in V gaute 
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between experiment and simulation below the kink allows for an estimate of the surface 
charge density. 



D. Holes 



As the electrochemical potential in the wire starts to approach the valence band, holes 
will appear. The creation and annihilation of holes are slow processes. The holes, like the 
surface states, will therefore not directly contribute to the capacitance at high frequencies 
but instead screen the electrons. The high effective mass of the holes imply a large density of 
states which results in a strong screening. If holes are created, practically no electrons will be 
removed any longer from the conduction band. The constant electron concentration causes 
the capacitance curve to flatten out at a finite value. This can be seen for the simulated 
doping Nd = 4 • 10 18 cm~ 3 in Fig. [3J The lower doping concentrations are too small for 
inversion to start before the last electrons have been removed, and the capacitance goes to 
zero in agreement with the experimental data. 

From Fig. H] we see that inversion has an effect for the thicker wires also for Np = 
2.2 • 10 18 cm -3 . This is in agreement with the measured capacitance which seems to flatten 
out at a finite value for negative bias rather than dropping to zero as for the thinner wires. 

In this context it is interesting to evaluate the doping required for the coexistence of 
holes and electrons. For zero temperature, this situation corresponds to the electrochemical 
potential entering the valence band at r = R nw before it leaves the conduction band at r = 0. 
Thus, a certain doping concentration is required if the inversion should have any effect on 
the capacitance. This doping concentration can be classically estimated with help of Gauss 
law 



The required doping as a function of i? nw can also be estimated with the Poisson-Schrodinger 
solver. The results are displayed in Fig. [5] Due to the confinement energy, the prediction of 
the simulation is higher, than for the classical estimate. The difference is most significant 
for small R nw . For large values of i? nw the classical prediction is equal to, or even larger than 
the model prediction. The reason is that the full bandgap is used in Eq. (jSj). In contrast, due 
to the thermal distribution, holes are already occupied slightly before the electrochemical 
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FIG. 5: Doping density of InAs NWs required for the coexistence of holes and electrons 

potential enters the valence band. From Fig. [5] we conclude that at the estimated doping 
concentration Np = 2.2 • 10 18 cm -3 , inversion will not affect the capacitance for R QW < 30 
nm in accordance with the findings of Fig. |H 

IV. COMPARISON BETWEEN DIFFERENT MATERIALS 



We will now show that equating the capacitance of a NW with its geometrical limit can 
lead to large overestimates. In Fig. [6] the theoretically calculated capacitances of InAs-, 
InSb, GaAs- and Si-NW:s are plotted. The effective masses of these materials are 0.023, 
0.0145, 0.063 and 0.33 [14] respectively. (For Si the density-of-state mass per conduction 
band valley was used and the valley degeneracy was explicitly taken into account for.) They 
are all surrounded by Hf02 and have the same geometry, R nw = 25nm and -R ga te = 35nm, i.e 
they have the same geometrical limit. For comparison the results from a Poisson-Thomas- 
Fermi solver have been plotted as well. In the Thomas-Fermi approximation, the electron 
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concentration at each point is determined by the distance between the conduction band 
edge and Ep. Multiplying the Fermi-Dirac distribution with the three-dimensional density 
of states and integrating over the conduction band yields the electron concentration as 

e{ ) J exp((E - e4>{r))f{kT)) + 1 1 ) 

where the conduction band edge is given by the electrostatic potential times the electron 
charge —e<p{r). Again we use Ep = as a reference. 

Due to the complex band structure of Si, only the Poisson-Thomas-Fermi results have 
been calculated for this material. The main difference between the two solvers is that 
classically the electron concentration does not have to be zero at the semiconductor-oxide 
interface, resulting in a higher capacitance. It is evident that the two models give the same 
results for a repulsive gate voltage. In this case the electron concentration approaches zero 
at i? nw as the potential surpasses the electrochemical potential, see also 121 ]. 

There is a reason why the geometrical capacitance cannot be reached even classically. 
Let us compare with a metallic coaxial cable. In this case the geometrical limit is reached 
as all charge is added at the surface. This requires an infinite density of states, which can 
not be found in a semiconductor, and thus the geometrical capacitance will not be reached. 
Materials with a high effective mass, such as Si, have a large density of states and will 
therefore come closer to the geometrical limit. 

Frequently, it is assumed that the capacitance is constant for large positive gate voltages. 
Both experimental and theoretical results point out that this is not the case, the capacitance 
is slowly increasing with bias. As the gate becomes more attractive, electrons will be filled 
higher up in the conduction band. Since the 3D density of states goes as E 1 / 2 , the density 
will increase as the bias is raised, resulting in a larger capacitance. 



V. CONCLUSIONS 



It has been shown how information about the doping and surface state concentrations 
can be obtained for semiconductor nanowires by comparing experimental and simulated 
capacitance data. The above procedure suggests how the contribution from these two factors 
can be separated. Analysis shows that the doping density Np> = 2.2 • 10 18 /cm 3 gives good 
agreement with the experimental data. The surface state charge at full depletion is estimated 
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FIG. 6: Comparison between theoretical capacitance and geometrical limit = 2.48 • 10 18 F/nm, 
for i? nw =25 nm i2 ga te =35 nm No = 2 • 10 18 /cm 3 , plotted against the electrostatic potential at 

to be around 10 13 e per cm 2 . Finally, the actual capacitance can be significantly smaller than 
its geometrical limit even in accumulation mode. The difference is of particular importance 
for nanowires with low effective mass. 
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